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Abstract. In terms of the stochastic process of quantum-mechanical version of Markov chain 
Monte Carlo method (the MCMC), we analytically derive macroscopically deterministic flow 
equations of order parameters such as spontaneous magnetization in infinite-range (d(— oo)- 
dimensional) quantum spin systems. By means of the Trotter decomposition, we consider the 
transition probability of Glauber-type dynamics of microscopic states for the corresponding 
(d + l)-dimensional classical system. Under the static approximation, differential equations 
with respect to macroscopic order parameters are explicitly obtained from the master equation 
that describes the microscopic-law. In the steady state, we show that the equations are identical 
to the saddle point equations for the equilibrium state of the same system. The equation for 
the dynamical Ising model is recovered in the classical limit. We also check the validity of 
the static approximation by making use of computer simulations for finite size systems and 
discuss several possible extensions of our approach to disordered spin systems for statistical- 
mechanical informatics. Especially, we shall use our procedure to evaluate the decoding process 
of Bayesian image restoration. With the assistance of the concept of dynamical replica theory 
(the DRT), we derive the zero-temperature flow equation of image restoration measure showing 
some 'non-monotonic' behaviour in its time evolution. 



1. Introduction 

In Bayesian statistics and statistical inference within the framework, to evaluate the expectation 
of microscopic labels (such as 'pixels' or 'bits') over the posterior distribution of these microscopic 
states or to calculate the marginal probability is very important procedure to solve various 
problems in the research field of 'massive' probabilistic information processing pQ. For most 
of the cases, it is very hard for us to carry out the task except for several limited cases of 
solvable probabilistic models. Recently, to overcome the difficulty, deterministic approaches 
based on mean-field approximation including the belief-propagation succeeded in dealing with 
the problems efficiently [21 [3] . 

On the other hand, in a lot of research fields, the Markov chain Monte Carlo method (the 
MCMC) has been widely used to calculate these expectations or to make marginal distributions 
by sampling the important states that contribute effectively to the macroscopic quantities 
[3]. Generally speaking, the MCMC takes a long time to wait until the Markovian stochastic 
process starts to generate the microscopic states from well-approximated posterior distribution. 
Especially, for some classes of probabilistic models which are categorized in the so-called spin 
glasses El \7\, the time consuming is sometimes very serious to make attempt to proceed 



the desired information processing within a realistic time. However, even for such cases, various 
improvements based on several important concepts have been proposed and succeeded in carrying 
out the numerical calculations [HJ [9]. From the view point of statistical physics, transitions 
between microscopic states are controlled by a specific hyper-parameter, namely, 'temperature' 
of the system and by cooling the temperature slowly enough during the Markovian stochastic 
process, one can get the lowest energy states efficiently. This type of optimization tool based on 
'thermal fluctuation' is referred to as simulated annealing 

Recently, the simulated annealing has been extended to the quantum-mechanical version. 
This new type of simulated annealing called as quantum annealing [T2J Q2J O E2 CES] is 
based on the adiabatic theorem of the quantum system that evolves according to Schrodinger 
equations. To use the quantum annealing, or more generally, to utilize the quantum fluctuation 
for combinatorial optimization problems or massive information processing dealing with a huge 
number of particles such as image restoration or error-correcting codes, the approach by solving 
the Schrodinger equations is apparently limited and we should look for different ways to simulate 
the quantum systems. 

As the most effective and efficient way, the quantum Monte Carlo method [T7] was established. 
The method is based on the following Suzuki- Trotter decomposition for non-commutative two 
operators A and B : 



Namely, the d-dimensional quantum system is mapped to the corresponding (d+ l)-dimensional 
classical spin systems. This approach is very powerful and a lot of researches have succeeded 
in exploring the quantum phases in strongly correlated quantum systems. However, when we 
simulate the quantum system at zero temperature in which quantum effect is essential, we 
encounter some technical difficulties although several sophisticated algorithms were proposed 
|18| . From the view point of information science, it is very informative for us to evaluate the 
process of information processing and if one seeks to utilize the quantum fluctuation to solve 
the problems, we should use the 'zero-temperature dynamics'. For this purpose, it seems that 
we need some tractable 'bench mark tests' to investigate the 'dynamical process' of information 
processing such as quantum annealing by using quantum fluctuation at zero temperature. 

In classical system, Coolen and Ruijgrok [TO] proposed a way to derive the differential 
equations with respect to order-parameters of the system from the microscopic master equations 
for the so-called Hopfield model [20J as an associative memory in which a finite number of 
patterns are embedded. The procedure was extended by Coolen and Sherrington [21], Coolen, 
Laughton and Sherrington [22] to more complicated spin systems categorized in the infinite- 
range (or mean-field) models including the Sherrington-Kirkpatrick spin glasses [23J. The so- 
called dynamical replica theory (DRT) was now established as a strong approach to investigate 
the dynamics in the classical disordered spin systems . On the considering the matter, it seems 
to be very important for us to extend their approach to quantum systems evolving stochastically 
according to the quantum Monte Carlo dynamics. 

In this paper, in terms of the stochastic process of quantum-mechanical version of Markov 
chain Monte Carlo method, we analytically derive macroscopically deterministic flow equations 
of order parameters such as spontaneous magnetization in infinite-range (d(= oo)-dimensional) 
quantum spin systems. By means of the Trotter decomposition, we consider the transition 
probability of Glauber-type dynamics of microscopic states for the corresponding (d + 1)- 
dimensional classical system. Under the static approximation, differential equations with respect 
to macroscopic order parameters are explicitly obtained from the master equation that describes 
the microscopic-law. In the steady state, we show that the equations are identical to the saddle 
point equations (equations of states) for the equilibrium state of the same system. We easily 





find that the equation for the dynamical Ising model is recovered in the classical limit. We also 
check the validity of the static approximation by computer simulations for finite size systems 
and discuss several possible extensions of our approach to disordered spin systems for statistical- 
mechanical informatics [23]. Especially, we shall use our procedure to evaluate the decoding 
process of Bayesian image restoration [23 E3 EH 1221 EH]. With the assistance of the concept 
of dynamical replica theory (the DRT) [21\ [22] , we derive the zero-temperature flow equation 
of image restoration measure, namely, overlap function between a given original image and the 
degraded one, showing some 'non-monotonic' behaviour in its time evolution. 

This paper is organized as follows. In Section 2, we explain our formulation to describe the 
deterministic flow equation from master equation for a simplest quantum spin system. In Section 
3, we apply our approach to evaluate the decoding process of Bayesian image restoration. Some 
'non-monotonic' behaviour in time evolution of image restoration measure is observed by the 
analysis of flow equations. The last section contains some remarks. 



2. The model system and formulation 

In this section, we derive the differential equations with respect to several order parameters for 
a simplest quantum spin system, namely, a class of the infinite range transverse Ising model 
|30t I3T] described by the following Hamiltonian: 
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where af and af denote the Pauli matrices given by 




It should be noticed that from the view point of Bayesian statistics, the above Hamiltonian 
corresponds to the logarithm of the posterior distribution. For various choices of parameters 
{Jij},h and {r}, one can model the problem of information processing appropriately. For 
instance, for the choice of Jjj = J, h ^ 0, the above Hamiltonian describes image restoration 
PU1 12"51 |2T)] from a given set of degraded images r = (ti, • • • , tjv) by means of Markov random 
fields. For Jy ^ 0, h = 0, it corresponds to the logarithm of the posterior of the Sourlas 
codes [32l [33] sending the parity check of two-body interactions through the binary 

symmetric channel (BSC). For the Sherrington-Kirkpatrick model [23j, we may choose Jij 
obeying the Gaussian with Jo mean and J 2 variance. Especially, for the choice of the so- 
called Hebb rule Jy = Y^u=i ^i^j' ^ = 0, it becomes energy function of the Hopfield model [20] 
in which extensive number of patterns (loading rate a) are embedded. 

In this section, we focus on the simplest case of Jij = J > 0, h = ferromagnetic transverse 
Ising model |34[ 135] . Let us start our argument from the effective Hamiltonian which is 
decomposed by the Suzuki- Trotter formula ([1]). 
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where k means the Trotter index and M denotes the number of the Trotter slices. We also 
defined the parameter B as B = (1/2) log coth(/3r/M) and a microscopic spin state on the k-th 
Trotter slice by cr fe = (oi(&), (T2(k), • • • , o"at(A;)), a- L {k) £ {+1, — 1}. 



2.1. The Glauber dynamics and its transition probability 

In the expression of the effective Hamiltonian ([3]), f3(fti(cr k , o~i{k +1)) is a local field on the cite 
i in the k-th. Trotter slice, which is explicitly given by 

j 

It should be noticed that the effective Hamiltonian §3§ is a classical system under the Trotter 
decomposition. Therefore, the quantum Monte Carlo dynamics should be described by the 
Glauber-type stochastic update rule whose transition probability is explicitly written by 

Wi(a k ) = - [1 - ai{k) tanh(/3(^j(tr fe : a(k + 1)))] . 



Namely, this means that 
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Obviously, in the classical limit, namely, when goes to infinity as T — > and (pi{cr k : <7j(fc+l)) ~ 
Bo~i(k+l), o~i(k) is identical to <7j(A;+l) with probability 1. Thus, for small T, we have cr k = tT^+i 
with a relatively high probability. In this sense, we can regard the term Bo~i(k+l) as an 'external 
field' from the nearest Trotter slice k + 1. Here we use the periodic boundary condition for the 
Trotter direction, that is, cry = cr^j + \. 

2.2. The master equation 

Then, the master equation for the probability of the microscopic states on the k-th. Trotter slice 
Pt{cr k ) is written by 
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where pt(cr k ) denotes a probability that the system in the A;-th Trotter slice is in a microscopic 
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state cr k at time t. We also defined an operator F£ to flip a single spin on the cite i in the k-th 
Trotter slice as <7j(fc) — > —o~i(k). 

2.3. From master equation to deterministic flow of order parameter 

We next introduce the spontaneous magnetization in the fc-th Trotter slice m k = iV _1 ^ o~i(k) 
as a relevant order parameter. The probability that the system is described by the magnetization 
rrik at time t is given in terms of the probability pt(cr k ) for a given realization of the microscopic 
state as Pt(mk) = Yler k Pt( a 'k)fi{™'k — mkip'kj)- Taking the derivative of this equation with 
respect to t and substituting ([5|) into the result, we obtain 

dPt ^^ = g~ ^Ptjfk) ^ ~4^- [1 ~ Qi(fc) tanh(/3<fo(<T fc : cn{k + 1)))] 5(m k - m k (a k )) 
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by neglecting 0(N 1 ) term. To continue the derivation of deterministic flow equations of order 
parameters, we use the following assumption: 

J im T7 y>nh[/3<Mo- fc :<ri(k + l))] = (tanh[/30(A;)]) Wfc) (7) 
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where we defined effective single site local field /3(f>(k) = ((3J/M)mk + Ba{k + 1) and the average 
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Namely, here we shall assume that in the thermodynamical limit N — > oo, the physical 
quantity in a particular choice (realization) of the Trotter slice, say, the fe-th Trotter slice 
iV -1 tanh[/3^>j(<Tfc : Oi(k + 1))] is identical to its own average over all possible paths in the 
imaginary-time axis for both 'past' {<r(l), <r(2), • • • ,cr(k— 1)} and 'future' {<r(A; + l), • • • ,cr(oo)}. 
The weight of each path is proportional to 
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In order to grasp the physical meaning of the quantity ((7|), it might be helpful for us to notice 
that tanh[/30(fc)] can be rewritten as 
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Thus, the quantity t&nh[f3<ft(k)] means a conditional expectation over the effective single spin 
u{k) (in the thermodynamic limit N — > oo) for a given 'past' {<r(l), • • • , a(k—l)} and it becomes 
a function of a(k + 1). After averaging over the 'past' with the weight P(a(l), ■ ■ ■ ,cr(k)) and 
'future' with the weight P(a(k + !),-•• , cr(oo)), we obtain 
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where we defined the summation of all possible paths by tr{ CT } (•••) = Eo-(i)=±i ' ' ' Eo-(M)=±i ('"')• 
We should notice that in the classical limit B — > oo (r — > 0), the weight of the path in the 
imaginary-time axis is dominated by <r(l) = • • • = a(M), mi = ■ ■ ■ = m as we always see it 
in the path integral approach of quantum mechanics in the limit of the zero Planck constant as 
h -¥ [361 EI] , namely, P(er(l), • • • , er(M)) = UiL <K°"(0 - o) and we have 
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Substituting the result into ([6]) and carrying out the integral with respect to rrik after multiplying 
the mfc, we immediately obtain the spontaneous magnetization flow for the dynamical Ising 



model; dm/dt = — m + tanh(/3 Jm) . Near the critical point, it behaves as dm/dt = — (1 — 
(3J)m — (/3 Jm) 3 /3 + 0(m 5 ). From this equation, we easily find that spontaneous magnetization 
shows well-known time-dependent behaviour m(t) ~ m(0) Q~ t K l ~P J ) around the critical point 
(3 ~ (3 C = J™ 1 , and at the critical point, the relaxation time diverges as (1 - /3 C J) _1 resulting m 
m(t) ~ t -1 / 2 (critical slowing down with dynamical exponent v = 1/2). 

Therefore quantum fluctuation comes from the finite B (T > 0). For the quantum case, the 
equation ^ leads to 

= J-{m k P t (m k )} - J-{P t (m k )(a(k)) path } (9) 
dt orrik onik 

In order to obtain the deterministic equation of order parameter, we should use the static 
approximation m k = m V(/c). Under this assumption and using the inverse procedure of the 
Suzuki- Trotter decomposition ([T|) : 
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we immediately have {a(k)) path = lim M ^oo(M 1 J2 k cr(k)) pat h as 
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and equation ([9]) leads to 
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Finally, substituting the form Pt(m) = S(m — m(t)) and making the integral by part with respect 
to m after multiplying itself m, we obtain the following deterministic equation. 

dm Jm 
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It is easy to see that the steady state dm/dt = is nothing but the equilibrium state described 
by the equation of state m = Jm{^J (Jm) 2 + T 2 } -1 tanh/3y (Jm) 2 + T 2 . In Figure QJleft), we 
plot the typical behaviour of zero-temperature dynamics (equation (|10p with /3 = oo) far from 
the critical point T c = J = 1 of quantum phase transition. We easily find that the dynamics 
exponentially converges to the steady state. The right panel denotes the zero-temperature 
dynamics at the critical point. The inset shows the log-log plot of m(t) indicating that the 
dynamical exponent in the critical slowing down is v = 1/2. This fact is directly confirmed from 
equation (jlpp with j3 = oo near the critical point dm/dt ~ —(1 — JT~ 1 )m — (Jm) 3 /2r 3 + 0(m 5 ), 
namely, m(t) behaves around the critical point r ~ T c = J as m(t) = m(0) e - *^ 1 "^ 1 " ' . At 
the critical point, the relaxation time diverges as tt = (1 — Jr" 1 )" 1 resulting in the critical 

slowing down as m(t) ~ e~ t ' TT — > m(t) ~ t~ v , v = 1/2. Of course, the exponent is the same 
as that of the 'mean-field model' universality class. 



Figure 1. Typical behaviour of zero-temperature dynamics described by (|10|) with ft — oo far from the critical 
point F c = J = 1 of quantum phase transition (left). The right panel denotes the zero-temperature dynamics at 
the critical point. The inset shows the log-log plot of m(t) indicating that the dynamical exponent in the critical 
slowing down is v — 1/2. 

2-4- On the validity of static approximation 

Without the static approximation, the following deterministic flow equations for each Trotter 
slice is obtained by substituting the form P t (mk) = S(mk — rrik(i)) and using the same way as 
deriving (fT0|) as follows. 

dm k tr W (j(/c) exp[^ ^fi x m^(/) + J B^fi 1 CT(/)fj(/ + 1)] 

— - — = — m k + hm — — Tn 57 7-7 (11) 

dt M— kxi tr w «p[gE£i^(0 + ^E£i^Mi + l)] 

Obviously, the equation ()lip is symmetric for the choice of k as long as we use the periodic 
boundary condition eri = (tm+i- This might be a justification to assume that the static 
approximation is correct at least for the present pure Ising system. 

To confirm this argument, we carry out computer simulation for finite size system having 
N = 400 spins. We observe the time evolving process of the histogram P(mfc) which is calculated 
from the M = N = 400 copies of the Trotter slices. We show the result in Figure [2j In this 
simulation, we chose the initial configuration in each Trotter slice randomly (we set each spin 
variable cr^ (fc) to +1 with a fixed probability p) and choose the inverse temperature (3 = 2 for 
r = 0.5 and T = 0.6. The time unit (the duration) of the update of P(mfe) is chosen as 1 Monte 
Carlo step (MCS). From both panels in Figure [2j we find that at the beginning, the P(m&) is 
distributed due to the random set-up of the initial configuration, however, the fluctuation rapidly 
(eventually) shrinks leading up to the delta function around MCS ~ 100. After that, the P(mfc) 
evolves as a delta function with the peak located at the value of spontaneous magnetization 
which is explicitly indicated in the inset of each panel. It should be noted that we evaluated 
the value of order parameter at the time point in the Runge-Kutta method. Thus, the duration 
between the points to be evaluated is not the MCS but the Runge-Kutta step. Of course, some 
statistical errors for the finite system should be taken into account, however, the limited result 
here seems to support the validity of the static approximation even in the dynamical process. 

3. Disordered systems: An application for Statistical-Mechanical Informatics 

It is easy for us to extend the above formulation to some class of disordered spin systems, that is 
to say, the infinite-range random field Ising model which is often used to check the performance 
of image restoration analytically [26 . 128t [29] as a bench mark test. 




Figure 2. Time evolution of the distribution P(mk) calculated for finite size system with N = M = 400. We 
choose the inverse temperature /3 = 2 for T = 0.5 (left) and F = 0.6 (right). The inset in each panel denotes the 
deterministic flows of spontaneous magnetization calculated by ffTTTf) for corresponding parameter sets. 



Here we consider a given original image £ = (£i, • • ■ , £/v) which is generated from the infinite- 
range ferromagnetic Ising model whose Gibbs measure (distribution for effective single spin) is 
described by -P(£i) = e /3sm °^/{2 cosh(/3 s mo)}, where mo denotes spontaneous magnetization at 
temperature f5j . A snapshot £ from the distribution is degraded by additive white Gaussian 
noise (AWGN) with mean ao£i and variance a 2 , namely, each pixel in the degraded image 
r = (ti, • • • , ttv) is obtained by r* = ao£i + axi with X{ ~ A/"(0, 1). From the Bayesian inference 
view point, we assume that the posterior P(<t\t) (here we define a as an estimate of the 
original image £) might be proportional to the logarithm of the effective Hamiltonian ([2]) with 
Jij = J V(i,j) and h ^ 0. The first and the second terms appearing in the right hand side of 
([2]) correspond to the prior distribution and the likelihood function, respectively. Whereas, the 
third term is introduced to utilize quantum fluctuation to construct the Bayes estimate for each 
pixel ('majority-vote decision' on each pixel), namely, sgn((af}). 

In this section, we attempt to describe the recovering process of original image through the 
deterministic flows of several relevant order-parameters and image restoration measure, namely, 
the overlap function M. = N^ 1 ^ £i s g n (( cr J z ))- 

For the above set-up of the problem, the local field on the site i in the fe-th Trotter slice now 
leads to 

<j)i{cr h : ai(k + l),r) = -j^^(Tj(k) + -^n(Ti(k)+Bai(k + l). 

3 

As relevant order parameters, we choose Hk = TV -1 TjCJj(fc) and magnetization m k . Then, we 
derive the differential equation with respect to Pt(m k , fj, k ) = ]C<x fc Pt( <T fc)<5( m fc — r71 -fc(""fc))^(/ u fc — 
fJ'ki^k)) as follows. 
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By assuming the self-averaging properties on the following physical quantities over both all 
possible paths in the imaginary-time axis and input data; original images and degrading 
processes (a particular realization of the quantity is identical to the average value and its 
deviation from the average eventually vanishes in the limit N — > oo), we have 
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where we defined the two different kinds of the averages by 
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Under the static approximation, we obtain 
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with Sm ,a °\i, x) = Jm + haoC + /iax and Dx = dx e~ x2 l 2 / \/2tt. Using the same way as the pure 
Ising system discussed in the previous section, we finally obtain the deterministic flow equations 
of the order-parameters m and n as follows. 
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For the solution of the above deterministic flows (m, fx) at time t, the overlap between the original 
image and degraded image is measured by 
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where (m, fi) is a solution of the following coupled equations 
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for a given point on the trajectory (m, fj,) at time t. To obtain the overlap function (|16p 
and (fT7|) (fT8]l . we used the concept of dynamical replica theory (the DRT) [2TJ 22], namely, 
'equipartitioning' and 'self-averaging' of the A4 during the evolution in time. As the derivation 
is a bit complicated, we shall show the detail in |Appcndix A| 

We solve the equations (|14 p(jl5p with (|16j) - ()18j) numerically and show the results in Figure 
El We choose the set of the parameters for the original image as /3.T 1 = 0.9 (too = 0.523) and 
ao = a = 1 which means that the corresponding optimal hyper-parameters are h = ao/a 2 = 1 
and J -1 = /3 . We consider the zero-temperature restoration dynamics in which the fluctuation 
to make the Bayesian estimate is only quantum-mechanical one (it is controlled by the amplitude 
of quantum- mechanical tunneling T). In the left panel, the deterministic trajectories in the space 
(to, y) are plotted for V = 0.6. The state of system in which the image restoration is successfully 
achieved is in ferromagnetic phase. Thus, order-parameters m and converge to the fixed point 
exponentially (there is no critical slowing down in this model system). In the right panel, we 
show the time evolution of the image restoration measure A4 for several values of T. From 
this panel, we find that some 'non-monotonic' behaviour is observed at the initial stage of the 
dynamics when we fail to set the amplitude to its optimal value (r ~ 0.6). Similar behaviour 
was reported in the Bayesian image restoration via thermal (classical) fluctuation [38J. 
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r = o.i 
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Figure 3. Trajectories in the phase space (m, /i) are plotted for V — 0.6 (left). The parameters are chosen as 
PJ 1 — 0.9 (mo = 0.523), ao = a = 1 and the corresponding hyper-parameters as h — ao/a 2 — 1 and J -1 = 
The right panel shows the evolution of image restoration measure, overlap function M in time for several values 
of the amplitude of quantum-mechanical tunneling F. 

In the Bayesian framework, it is desired for us to obtain the estimate of each pixel and 
hyper-parameters simultaneously. In such case, we use the so-called EM algorithm based on 
the maximization of marginal likelihood criteria [39j. The quantum-mechanical extension and 
the formulation presented here is applicable to the simultaneous estimation for both micro and 
macro parameters. 

4. Concluding remarks 

In this paper, for a simplest quantum spin systems, we showed a formulation to describe 
the macroscopically deterministic flows of order parameters from the master equation whose 
transition probability is given by the Glauber-type. Under the static approximation, differential 
equations with respect to macroscopic order parameters were explicitly obtained from the master 
equation describing the microscopic-law. In the steady state, we found that the equations are 
identical to the saddle point equations for the equilibrium state of the same system. We also 
checked the validity of the static approximation by computer simulations and found that the 



result supports the validity of the approximation. Several possible extensions of our approach 
to disordered spin systems for statistical-mechanical informatics was discussed. Especially, we 
used our procedure to evaluate the decoding process of Bayesian image restoration. With the 
assistance of the concept of dynamical replica theory (the DRT), we derived the zero-temperature 
flow equation of image restoration measure showing some 'non-monotonic' behaviour in its time 
evolution. Of course, by using the present approach, one can evaluate the 'inhomogeneous' 
Markovian stochastic process of quantum Monte Carlo method (in which amplitude T is time- 
dependent) such as quantum annealing. In the next step of the present study, we are planning 
to extend this formulation to the probabilistic information processing described by spin glasses 
such as quantum Hopfield model |4U[ 141] including a peculiar type of antiferromagnet |42j . 



Appendix A. Derivation of overlap function: A dynamical approach 

In this appendix, we show the derivation of overlap function for image restoration problem 
()16p -(|18 p by using a useful concept of dynamical replica theory (the DRT) |2H I22j. 

The problem here comes from the fact that the quantity M. = N^ 1 Ylii^i s S n (( a i{k))path)]data 
containing a single site expectation ((Ji(k)) pat f l is not observable in computer simulations of a 
single system. Therefore, we must introduce infinite number of 'virtual copies' ('real replicas' 
or 'ensembles' ) to evaluate the overlap function by using the law of large number. Namely, we 
evaluate the 'real-time' dependence of the overlap function as follows. 



M(t) 



copy- 



data 



M(m k (t),n k (t)) 



where p = 1, • • • , L denote the copies of the original system in which pt(a p {k)) is the same. From 
the law of large number, we can expect lim/^oo L^ 1 X^p=i a lt(^) = ( a i(k))copy To calculate the 
expectation (• • • ) C opy, we assume the equipartitioning in the (m^, /ifc)-subsells. Namely, explicit 
time-dependence of the expectation through pt{a p {k)) is now removed within the subsells as 



(...) 



E {(T P( fc ) } np=iPt(g p (fc))(- • • )5(m k - m k (*P(k)))5{n k - ^(k))) 

£ { <TP( fc)} Ilp=i Pt(trP(k))S(m k - m k (<rP(k)))6(n k - fj, k ((rP(k))) 
E{ CT p(fc)} Up=i(- ■ ■ )${ m k ~ m k {aP{k)))5{n k - n k (crP(k))) 



oo as 



E{o-p(fc)} Up=i 5(m k - ra k {aP(k)))5{^ k - ^ k {aP(k))) 
Then, assuming the self- averaging on the we have M.i,{m k ^\x k ) in the limit of iV 
M L {m k ,n k ) 

, Eo-i(fc),- ,<r L {k) FlJ=i s ( m k ~ m k {aP{k)))5{^ k - n k (crP (fc)))&sgn (± £^ =1 a p (k) 

lim 

7V-s>oo 



E 



<T 1 (k),— ,cr^{k) llp=i S(m k - m k (aP(k)))5(fi k - fi k (aP(k))) 



lim lim 

n->0 AT->oo 



( L \ n L 

iE^w nn*( 
p=\ ) a=\p=\ 



m k - m k (*P a (k)))5(n k - v k (* pa (k))) 



where we used the fact 
"tr {s }PV" 



tr {s} P 



data 



lim 

n— >0 



(tv {s} vmH S }'p) n 



data 



lim 

n-s-0 



a=l 



data 



and introduced the replica index a = 1, • ■ ■ , n to carry out the average [• • • }data by standard 
replica trick. It should be noted that a = (eri (&),•■■ , <t/v(&)) and {• • • } denotes the set of 
spin variables for all possible combinations of copies and replicas (p, a). By using the integral 
representation for the delta-function, we obtain 



MiK.ft) = lim lim TT e -^ pa (m k fn^ + ^n 

J - QO pa 

x £ 6sgn(iX;<(A;)]e-^p^rE^rW 



data 



Now, under the static approximation, the quantum fluctuation appears through the quantities 
= m and fj,^ = fx which obeys (|14j) (|15|) at any time t. Therefore, from now on, we cancel the 
/c-dependence. By simple transformation of the variables as —ip pa h-> p pa ,—im pa \— > m pa and 
carrying out the data average [• • • }data over the effective single site distribution (fT3|) . one obtains 



M L (m,p) = lim lim f°° TT e ^E P J^r+^P 



2 cosh(/3 s mo) 



6sgn (1 Ej=i *f ) exp[E pQ m^a^ + « + ax) £ pa A pQ ^ Q ] 
£ ct p* expE pQ mP« + K + ax) £ pQ A pf V^] 

/oo 
drhdp expL/VvHrn, 
-oo 

E. pa ftsgn (l gU *f ) exp[E pa + « + ax) £ pa ^ pa ] 

E ct p« exp[E pQ mP« + « + ax) £ pQ ^ 



with 

V(rh,p) = -^2(mm pa + fifi pa ) + log 



pa 



/ ni 2cosh K° + (aoC + ax)AHI 

2cosh(/3 s m ) ./-oo _ 



We should notice that in the above expression of the overlap function A^z,(m,//), the condition 
£i = sgn(L _1 E P =i °>i(A;)) leading to 'perfect image restoration', namely, limf^ 00 .A/fL = 1 
immediately gives lim n _>.o lhnjv->o I^°oo drhdp exp[N^(rh, p)] = 1. Therefore, it is easy to find 
that the part lim^^o linijV-^o J^L, drhdp exp[A r ^'(m, p)] appearing in (|A.1|> is a normalization 
factor. Thus, the overlap function derived by the above dynamical approach now leads to 

E„r« a sgn (i E P= i of) exp[E pa m"*^ + « + ax) £ pQ P"*^ 
M.L{m,p) = 



E ct p« exp[^ pa m p *a pa + « + ax) £ pQ P pa o pa ] 



where rh or p should be chosen as a saddle point of the function *&(m, p). Assuming the replica 
symmetric and the copy symmetric solution rh pa = rh,fi pa = p V(p, a), we obtain the function 
to be optimized at the saddle point. 



^rs = hni hm -! — = —mrh — uu-\ ^— —- / Dx cosh[m + (ao£ + ax)p] 

n-s-0 JV,L->oo IlL 2cOSh(p s ?72-o) J —oo 



Obviously we find that the saddle point is obtained by the equations (fT7|) and (fT8|) . Then, the 
overlap function is evaluated as 

MiXm.fi) 



I 00 Dx V sen f 1 V/ "TT * im+[a °' +axmul 
2 cosh(/3 s m ) Loo , I L 1 1 / 1 1 2 cosh [™ + (°o£ + ax)p] 

/oo / 1 i 

Dx ^ sgn(z)5 z - - ^ af 1 



00 ^ £ gfemo^ 

ciz ~ 

.oo 2cosh(/3 s m ) 7_ 



e 

x 



with 



n- — — -. — —7 = / dzsgn(z)VL(z) 
x 2cosh[m + (a £ + ax)p] J ^ 6 v ' w 



2cosh(/3 s m ) 7_oo , ./-oo 2 ^ ^ 2cosh[m + (a £ + ax)p] 

a 1 ,— > (T U P— 1 

/°° dy^ Eg£e^ ^ f cosh[m + (a g + ax)A - i^ 1 ?/] j L 
7_oo 2?r 2cosh(/3 s m ) 7_oo I cosh[m + (a £ + ax)/}] J 

= / T-j V * \(o \ / ^{cos(yL- 1 )-isin(yL- 1 )tanh[m + (aoe + ax)A]} 

J-oo 2 ^ 2cosh(/^ s m ) J_oo 

Thus, we have 

P(z) = Um P L (z) = \ , n Dx e iy(z-t^h[rh+(a^+ax)iA) 

L->oo 2cosh(p s mo) 7_oo J— 00 

2vr 



/oo 
Dx 5(z — tanh[m + (ao£ + ax)jj\). 
-00 



2 cosh(/3 s mo) 

Substituting the result into A4l and taking the limit of L — )■ 00, we finally obtain 

/OO /'OO 
(izsgn(z) lim Vl{z) = / cfe sgn(z)V(z) 
-00 ./-oo 



2 cosh (/3 S mo) 



Dxsgn[m + (ao£ + ax)/}] 



where we used the fact sgn(tanh(x)) = sgn(x). This result is nothing but equation (|16p . As 
parameters m and (1 are related to the order-parameters m and ^ in equations (|17j) (|18j) , overlap 
function M. is influenced by quantum fluctuation through (fT7|) fp~8|) and the solution of equations 
(EH USD. 
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